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Abstract: The Taylor expansion method has been used to investigate the scale depen- 
dence of the power spectrum of the curvature perturbation. In the present study, an alter- 
native numerical method is used to clarify the k dependence. Although there is thought 
to be no large difference between these two methods, some differences arise among various 
inflation models. For example, at k = 1 Mpc, there is a 1.4 % difference in the power spec- 
trum, and with respect to the angular power spectrum, the difference of the value of y 2 
nearly 10 occur in new inflation. However, in hybrid inflation, these differences do not oc- 
cur. The time dependence of the inflationary and cosmological parameters is investigated, 
and differences among inflation models are clarified. 
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1. Introduction 



The curvature perturbation produced during the inflation epoch is an important quantity in 
cosmology. The scale dependence of the power spectrum of the curvature perturbation has 
been investigated in detail using the Taylor expansion around a pivot scale [1]. The scale 
dependence of the power spectrum can be written in terms of the slow-roll parameters 
and can be calculated in terms of higher orders of the slow-roll parameters [2]. In the 
present paper, a numerical method is proposed as an alternative method by which to 
investigate the k dependence of the power spectrum. The proposed method is described 
in detail herein. In the case of a numerical method, the time dependence of the inflaton 
must be known in order to investigate the k dependence of the power spectrum. In a 
previous study, we investigated the effect of the length of inflation [3]. However, in the 
present study, we assume that the length of inflation is very long, and so do not consider 
the contribution of the length of inflation. The time-dependent behavior of the inflaton 
during the last 100 e-folds in inflation is investigated using simple inflation models, because 
the physically interesting period of inflation is that in the last 60 e-folds. Although a 
correct inflationary potential based on the superstring or supergravity theory has yet to 
be established, the potential term in this latter period of inflation can be expressed in a 
simple form, such as the new inflation model, the chaotic inflation model, or the hybrid 
inflation model. Here, these three inflation models are considered in the present study. 
Next, using the derived time-dependent behavior of the inflaton, the time dependence of 
the slow-roll parameters, the Hubble parameter, and the spectral index can be calculated. 
The difference in time-dependent behavior of such parameters among the inflation models 
can then be investigated. Next, using the time dependence of the slow-roll parameters 
and the Hubble parameter, the differential equation of the gauge potential is solved by 
a numerical method. Finally, the k dependence of the power spectrum of the curvature 
perturbation is derived using this solution. In order to investigate the difference between 
the familiar Taylor expansion and the newly proposed numerical method, the angular power 
spectrum and the value of \ 2 are calculated. 
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2. Cosmological and inflationary parameters 



The dependence of the cosmological and inflationary parameters on time is investigated 
using a slow-roll inflation model. Assuming a spatially flat universe and the energy density 
of the universe to be dominated by the inflaton field, the Einstein field equations can be 
written as 

and 

4> + 3Hj> + V'(<f>) = 0, (2.2) 

where 4> is the inflaton field, V{4>) is the inflaton potential, m is the Planck mass, and 
H = a/a with a being a scale factor. The variable ./V which is defined by dN = Hdt, is 
introduced to represent time with respect to the number of e-folds of inflation. Equation 
(2.1) can then be written using the relation <p = Hdcp/dN as 

3m 2 (l- 4ir/3m 2 (d^/dN) 2 )' 1 J 

The Hubble parameter is usually considered to be a constant, i.e., H 2 = 8n /3m 2 V((j)) . In 
the generalization, however, with (f> being time-dependent, the Hubble parameter is also 
time-dependent (from equation (2.3)). Using the variable N, equation (2.2) can be written 
as 

d 2 4> l_dtf_# M_ V'{4>) _ 

dN 2 + HdNdN +3 dN + ^I 2 ~~ { ' 

Equation (2.4) can be further written in terms of and the derivative of (ft with respect to 
N. In order to solve equation (2.4) numerically, three slow-roll inflation models are adopted: 
the new inflation model with the potential term given by V(cf>) = A 2 v 4 (l — 2((f>/v) p , (p = 
3,4), the chaotic inflation model with the potential term given by V(<j>) = M 4 /2((j)/m) a 
, (a = 2,4,6), and the hybrid model V(<j>) = a[(v 2 - a 2 ) 2 + m 2 /2(f) 2 + g 2 (f> 2 cr 4 } ~ a(v 4 + 
m 2 /2(f) 2 ). For example, in the case of chaotic inflation (a = 2) using an initial condition of 
N = no, (f>(no) = y 7 — n§/\[2nim and d<fi(no) / dN = —m/y/8ir(—no), the numerical solutions 
to equation (2.4) for no values ranging from -115 to -70 do not reveal any appreciable 
dependence on no, and, with the exception of the new inflation model (p = 3), this property 
can be shown for other inflation models. The behavior of <p for the case in which no = —100 
is shown in Figure 1 for three inflation models. The inflaton field slowly decreases in 
the cases of the chaotic inflation model and the hybrid model, but increases very slowly 
in the case of the new inflation model. The difference of the ^-dependence (time) of 
the Hubble parameter among the three inflation models is derived, and the behavior of 
H 2 (N)/H 2 (—60) is shown in Figure 2. There is no ^-dependence in the case of the new 
inflation model, but the Hubble parameter becomes considerably smaller as N approaches 
zero in the case of the other models: For example, the value of H 2 (N)/H 2 (— 60) is 1.33 for 
N = -80, 0.835 for N = -50, 0.669 for N = -40, and 0.504 for N = -30, in the case of 
chaotic inflation (a = 2). 
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The following parameters are used of slow-roll inflation [4] : 



4 ( £ + v)-' = £ ( *W)» p.) 



2 v 2 y 4vr v #(0) 

m 2 ff"(<ft) 
4V #(0) ' 



(2.6) 



m 4 H'WW'W 
^ 16vr 2 (#(</>)) 2 ' K ' 

Other slow-roll parameters (ey, J?y , £v) can written in terms of the slow-roll parameters 
e, 5, and £ to the first order in slow roll: e = ey, 5 = r/y — ey, and £ = £y — 3ey?iy + 3ey, 
where ey = m 2 /16tt(V /V) 2 , r/y = m 2 /87r(V7V), and £y = m 4 /64w 2 (V'V'" /V 2 ). The 
behavior of slow-roll parameters e, S, and £ with respect to the ^-dependence differs among 
the inflation models. 

The time dependence of the spectral index n s can also be calculated using these slow- 
roll parameters. The spectral index is given by 

n s = 1 + 25 - 4e. (2.8) 

The time dependence of n s is shown in Figure 3 for three inflation models. Figure 3 shows 
that, for each model, the decreasing behavior of the spectral indexes as well as their values 
and shapes are different. Taking the new inflation (p = 3) as an example, the value of n s 
decreases as N approaches zero: n s (— 80) = 0.950,n s (— 60) = 0.935,n s (— 50) = 0.921, and 
n s (-40) = 0.906 . 
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3. Scalar perturbations 



The dependence of the curvature perturbation on the parameter k is investigated assuming 
a spatially flat Friedman-Robertson- Walker (FRW) universe with a background spectrum 
given by metric perturbations. The line element for the background and perturbations is 
generally expressed as [5] 

ds 2 = a 2 (r){(l + 2A)dr 2 - 2d i Bdx i dr - [(1 - 2*)<% + 2d i d j E + h^dx^x^, (3.1) 

where r is the conformal time, the functions A, B, ^ , and E represent the scalar perturba- 
tions, and hij represents tensor perturbations. The density perturbations in terms of the 
intrinsic curvature perturbation of comoving hypersurfaces is given by TZ = — (H/4>)5(f), 
where cf> is the inflaton field, dcfi is the fluctuation of the inflaton field, and TZ is the curvature 
perturbation. Overdots represent derivatives with respect to time t, and primes represent 
derivatives with respect to the conformal time r. Introducing the gauge-invariant potential 
u = a(r)((5^ + ((j)/H)fy) allows the action for scalar perturbations to be written as [6] 

S= 1 -JdTdH(^) 2 -(Vuf + ^u 2 }, (3.2) 

where Z = a<j)/H and u = —Z7Z. The field u(t,x) is expressed using annihilation and 
creation operators as 

«(t,x) = J rf 3 M^(r)a k + <(r)a t _ k }e- kx , (3.3) 

and the field equation for Ufc(r) is derived as 

+ (k ~ = 0, (3.4) 

where the solution to satisfies the normalization condition u^du^/dr — u* k duk/dr = i. 
Using the slow-roll parameters, (d 2 Z/dr 2 )/Z is written exactly as 

^f-*"*V (3.5) 

where 

v = {l + e-h + e 2 -2e5+^+^-). (3.6) 
A partial integration allows r to be expanded as follows [7]. 

1 2e(e - 5) | f e(-25 2 + 5(9 - 4e)e - 6e 2 + 3e 3 - g(l - e)) dN (g ?) 



aH{l-e) aH{l-e) J (1 - e) 2 aH 

Here, using the approximate relation a 2 H 2 = (1 + 2e(e — 5)) 2 /r 2 (l — e) 2 , equation (3.4) 
can be rewritten as 

^C-^V. (3-8) 
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If it is assumed that the slow-roll parameters can be constant, the solution for equation 
(3.8) can be written using the Hankel function Hjp(—kr) as 



Ufc ( r ) = il e i( (1+ i/2),/ 2( _ jir) i/2 F M ( _ jtr)) (3 9) 

2y k 



where \i = \Jlv(\ + 2e(e — <5)) 2 /(l — e) 2 + 1/4. This solution is fixed such that as kr — > 
— oo, Ufe(r) approach plane waves. The numerical solution to equation (3.4) can then be 
considered. Equation (3.4) can be written in terms of N as 

d 2 u k 1 daH du k k 2 

dFP + ^H^NdN + W " >k = °- (3 0) 

In order to calculate equation (3.10) numerically, the present-day size perturbation k = 
0.002(1/Mpc) is assumed to exceed the Hubble radius in inflation at the time of N = —60. 
The scale factor is thus written as a = 0.002/ H Q e N+m , where H is the value of H at 
N = —60. Since the Hubble parameter and v have both been calculated as functions of 
N, equation (3.10) can be calculated numerically with an initial value of Uk derived from 
equation (3.9). Note that the desired solution of Uk is the solution at the point aH = k. 
Taking new inflation (p = 3) as an example, in Figure 4, the absolute value of Uk at aH = k 
derived by numerical calculation of equation (3.10) is shown as a function of k. The absolute 
value of Uk (equation (3.9)) obtained using the Bessel approximation, adopting the value 
of jjL as the numerically calculated value, is also shown. The value of changes slightly 
according to N: /i(-80) = 1.525, a*(-60) = 1.532, /x(-50) = 1.540, /i(-40) = 1.550 and 
^(—30) = 1.564. The results shown in Figure 4 do not suggest any appreciable difference 
between the numerical calculation and the Bessel approximations, indicating that the Bessel 
approximation, in which fj, exhibits k dependence, can be considered satisfactory. 
The power spectrum of the scalar perturbations Pr, is defined as follows [4]. 

<n k (T),n;(T)>=^-PnS 3 (k-i), (3.11) 

where, TZk{ T ) is the Fourier series of the curvature perturbation 1Z. The power spectrum 
Pti is then written as [4] 

where Z = a<j>/H. The power spectrum is estimated by three different methods: the 
numerical method, the Bessel approximation using equation (3.9), and the Taylor expansion 
with running (as a familiar method). Calculation of the power spectrum is estimated at 
k = aH. In the case of the Bessel approximation using the values of /x, d(f>/dN, and 
H derived from the numerical calculation (having k dependence), the power spectrum 
(equation (3.12)) can thus be written as 

_ k 2 (-fcrXfff^-fcr)) 2 , 
Pn ~ 8^ a^/dN) 2 lk=aH > (3 - 13) 

where the value of P-ji is obtained at k = aH (i.e. kr = —(1 + 2e(e — — e))). 

There is very little difference between the results of the numerical calculation and the 
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Bessel approximation (see Table 1), again demonstrating that the Bessel approximation is 
satisfactory. 

The power spectrum is commonly written as 

p * = ^ 1 O v - 1/2)I "^ 1 "'- (3 ' 14> 

The value of the power spectrum P-ji given by equation (3.14) is 4.0395 x 10 6 A 2 , assuming 
the value of fx at N = —60. Using the same parameter values, the numerical calculation 
by equation (3.12) and the Bessel approximation (equation (3.13)) both afford a value 
8.5488 x 10 6 A 2 . The above values are for the case of the new inflation (p = 3), but other 
inflation models exhibit the same behavior. There is thus a difference of approximately 
a factor of two between the present numerical and Bessel calculations and the familiar 
case (equation (3.14)). This discrepancy can be attributed to the difference in the Hankel 
function used to derive the expressions. The familiar expressions are derived from the 
asymptotic form (k/aH — > 0) of the Hankel function and are estimated at 4k/aH = 1 [4], 
whereas the present Bessel calculation uses the correct form. The difference between the 
asymptotic and correct forms of the Hankel function equates to a factor of approximately 
\/2 at k/aH = 1, giving rise to a change in the normalization of the power spectrum. 

Since the value of the spectral index differs from unity, as indicated by experimental 
results such as WMAP data [1,8], the power spectrum can be inferred to have k dependence. 
The power spectrum is usually expanded at a pivot scale ko, i.e. [1], 

P n (k) = i^(jfc )(A)n.(*o)-l+(l/2)aln(fc/*o) ( 3 . 15 ) 

ko 

where a = dn s /dlnk\k=k - Since the three inflation models have non-zero values of a, 
the running case must be considered. If this running term is ignored, with the exception 
of hybrid inflation, differences with respect to the power spectrum and the angular power 
spectrum occur between the cases with and without running. Taking the new inflation 
(p = 3) as an example, Table 1 lists the results for the three calculations for various 
values of k, and Figure 5 compares the numerically calculated power spectrum with that 
determined by the Taylor expansion with power spectrum running, where the overall factor 
Pn(ko) is fixed at 8.5488 x 10 6 A 2 . These results seem to indicate that the power spectrum 
given by the Taylor expansion exhibits behavior that differs from that obtained by the 
present numerical calculation in the case of large k. Here, the differences in k dependence 
of the power spectrum between the Taylor expansion and the proposed numerical method 
among the inflation models are as follows. A large difference, 1.3 % (p = 3) and 0.8 % 
(p = 4) occurs at k = 1 (1/Mpc) in the case of new inflation. In the case of chaotic 
inflation, the difference is 0.04 % in the 4> 2 case, 0.3 % in the </> 4 case, and 0.9 % in the </> 6 
case. The difference in the case of the hybrid model is slight. 

Finally, the angular power spectra are calculated using a modified CMBFAST code 
[9], and likelihood analysis is performed using the WMAP three- year data and the WMAP 
Likelihood Code [10]. Here, the value of \ 2 is calculated for these inflation models using 
the same values of the cosmological and inflationary parameters, but only the difference 
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of the value of x 2 is considered. For the new inflation model, the difference of the value 
of x 2 between the familiar Taylor expansion method and the proposed numerical method 
is nearly 10 for the case in which p = 3 and is approximately 3.8 for the case in which 
p = 4. For the chaotic inflation model, the difference of the value of x 2 between the familiar 
Taylor expansion method and the proposed numerical method is approximately zero at (ft 2 , 
approximately 1.4 at 4 , and approximately 2.7 at 6 . However, for the hybrid inflation 
model, the difference of the value of x 2 does not occur. The dependence of the initial value 
riQ in the numerical calculation in the range from -115 to -70 reveals no clear difference, 
except for new inflation (p = 3), where the difference of x 2 is approximately 2. 
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4. Discussion and summary 

We investigate the difference in k dependence of the power spectrum of the curvature per- 
turbation between the familiar Taylor expansion method and a numerical method, and 
time dependence of the cosmological and inflationary parameters using three slow-roll in- 
flation models. Differences in behavior are observed among the three inflationary models, 
i.e., the Hubble parameter changes only slightly with time (N) for the new inflation but 
decreases with time for the chaotic inflation model and the hybrid inflation model. The 
spectral index shows a similar decreasing behavior, but the value and shape are different. 
The calculation of the power spectrum using the Bessel approximation was shown to be 
consistent with the numerical calculations, and an approximately two-fold difference in 
normalization was shown to exist between the present methods and the usual methods (see 
equation (3.14)) due to the use of the asymptotic form of the Hankel function in the usual 
treatments. In the fc-dependent spectra, differences in the power spectrum and the angular 
power spectrum between the Taylor expansion method and the proposed numerical method 
occur among the inflation models. Specifically, in the new inflation, a large difference of 
Ax 2 — 10 occurs. In contrast, in chaotic inflation {(ffi and c/> 4 ), some difference of x 2 oc- 
curs, and in hybrid inflation and chaotic inflation {(f) 2 ), the difference of x 2 rarely occurs. 
As in the case of the new inflation model (p = 3), a number of unique and characteristic 
behaviors, as compared to other inflation models, are derived: 

1. The Hubble parameter changes only slightly in the range from N = —110 to N = —20. 

2. In the numerical calculation, there is a small no dependence of the power spectrum 
and angular power spectrum, but in other inflation models rarely occurs. 

3. The value of e (see equation (2.5)) is very small, and the value of |n| is larger than 
the value of e. 

4. The value of £, which is not very small, is larger than e and oscillates. Moreover, 
dependence of the initial condition no exists. 

Unfortunately, the WMAP three-year data was used in the calculation of x 2 ■ Since the 
difference of the value of x 2 is nearly 10, and the values of x 2 are larger than those for other 
inflation models in the case of the cosmological and inflationary parameters used herein, it 
is very interesting what a difference of the best cosmological and inflationary parameters 
between the Taylor expansion method and the proposed numerical method occurs using the 
WMAP five-year data. Moreover, we considered only the case of the scalar perturbation. 
In the future, we intend to investigate the effect of the tensor perturbation. 
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Table 1: Power spectrum values {P n x 10 6 A 2 ) and ratios of these values to those obtained through 
numerical calculation 



fc(l/Mpc) 


Numerical 


Bessel 


Taylor 


0.001 


8.9425 


8.9413 


8.9507 


(ratio) 


(1) 


(0.99987) 


(1.0009) 


0.002 


8.5488 


8.5488 


8.5488 


0.01 


7.7030 


7.7011 


7.6680 


(ratio) 


(1) 


(0.99974) 


(0.99544) 


0.1 


6.5924 


6.5907 


6.5305 


(ratio) 


(1) 


(0.99974) 


(0.99061) 


1.0 


5.6065 


5.6052 


5.5291 


(ratio) 


(1) 


(0.99976) 


(0.98619) 
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Figure 1: Variation in <p/m (m is Planck mass) with respect to N . The solid line repre- 
sents new inflation (p = 3), the dashed line represents chaotic inflation (a = 2), and the 
dash-dotted line represents hybrid inflation. 
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Figure 2: Variation in the Hubble parameter (H 2 (N)/H 2 (— 60)) with respect to N. The 
solid line represents new inflation (p = 3), the dashed line represents chaotic inflation 
(a = 2), and the dash-dotted line represents hybrid inflation. 
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Figure 3: Variation in the spectral index n s with respect to N. The solid line represents new 
inflation (p = 3), the dashed line represents chaotic inflation (a = 2), and the dash-dotted 
line represents hybrid inflation. 
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Figure 4: Variation in with respect to k. The dotted line represents the numerical calcu- 
lation, and solid line represents the Bessel approximation using the numerically calculated 
value of // in the new inflation. 
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Figure 5: Variation in Pji with respect to k. The dotted line represents the numerical calcu- 
lation, and the solid line represents the Taylor approximation with spectral index running 
for the same value of P-jz(ko) in the new inflation. 
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